The Zero Temperature Phase Diagram of the Kitaev Model 
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We show that the zero temperature phase diagram of the vortex free sector of the Kitaev model is 
in one to one correspondence with that of the classical dimer model on the same lattice. We find that 
the model generically has three distinct phases. On a honeycomb lattice with a 3 x 3 fundamental 
domain all three phases are accessible. As the couplings are varied there are two distinct transitions. 
The new transition is one to a gapped phase that opens up in the interior of the B phase. 
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Introduction. — Kitaev's observation [l[ that a spin 1/2 
system on a honeycomb lattice has a gapless phase with 
vortex excitations that obey non-Abelian statistics has 
stirred a lot of interest 0, S, H S, @, 0, B H- Yao 
and Kivelson || extended Kitaev's considerations to the 
Fisher (triangle-honeycomb) lattice and found that its 
zero-temperature ground state is a chiral spin liquid. We 
extend these considerations to the more general setting. 

Kitaev models divide into two classes: those defined on 
bipartite and non-bipartite lattices. A lattice is bipartite 
if its sites can be coloured black and white so that ad- 
jacent sites are always of the opposite colour. A Kitaev 
model is then bipartite if its Hamiltonian has no inter- 
action between bonds of the same colour. We establish 
that the vortex free, zero temperature, phases of Kitaev 
models are in one to one correspondence with the phases 
of classical dimer models. 

• Bipartite models have generically three phases, the 
third being a previously unnoticed gapped phase 
which we refer to as the C phase. This arises only 
when a sufficiently large fundamental domain is 
considered, e.g. a 3 x 3 hexagonal domain. 

As the couplings are varied the C phase opens up in 
the interior of the B phase and corresponds to giving a 
Dirac mass to the pair of degenerate massless Majorana 
Fcrmions of the B phase. Close to the transition when 
the mass gap is small the model is well approximated by 
a continuum massive Dirac field theory. 

• Non-bipartite couplings lift the degeneracy of the 
Majorana modes and the mass gap of one of the 
Majorana fields can go to zero to give a chiral spin 
liquid as discussed in 



After summarizing the Kitaev model and its formula- 
tion in terms of Majorana operators we establish that 
the quadratic form associated with the Hamiltonian in 
the vortex free phase is precisely the Kasteleyn matrix 
of a classical dimer model. This allows us to identify the 
phase diagrams of the two models and establish the ex- 
istence of a new gapped C phase in the bipartite Kitaev 
model and predict the mass gap. It also affords us an un- 
derstanding of non-bipartite models such as the extended 
Kitaev model of [f| where the couplings can be choosen to 
correspond to those of a classical Ising model in its dimer 
realization . At the Ising critical point the model has 
a low- lying massless spectrum that corresponds to a "rel- 
ativistic" chiral Fermion [ll| which becomes massive as 
the couplings are varied. 

The model. — The Kitaev spin model on the hexagonal 
lattice is given by the Hamiltonian 

n Y ■><<< ■ E E J *°M w 

a: — link y — link z — link 

and partition function 

Z = Tr{e~ m ). 



(2) 



For simplicity we only consider the case of open boundary 
conditions as our focus is on the phase diagram of such 
models. 

With open boundary conditions the model (p} can be 
rewritten, after a Jordan- Wigner transformation [J 0, 
5], in terms of Majorana operators Cj with the partition 
function © given by 



Z 



and the Hamiltonian 
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where {ci,Cj} = 2<5y. The matrix element Kij between 
site i and j is given by one of the couplings J x , J y or J z , 
up to an overall sign. This rewriting is valid for quite 
general couplings which can vary throughout the lattice. 
Yao and Kivelson [f| extended Kitaev's construction to 
the Fisher lattice (or triangle-honeycomb lattice). 

The most general Kitaev model then assigns positive 
couplings to all links and a Z2 gauge field Uij to each 
z-link between sites i and j. 

Unfortunately, for general temperatures the Kitaev 
model has not yet yielded an exact solution. However, 
the zero temperature ground state is known, thanks to 
Lieb's theorem [T3|, to be the vortex free state. 

Phase diagram equivalence. — In the vortex free case 
the Uij become w|j d , and is the Kasteleyn matrix 
of a classical dimer model. In this context uf^ is a 
Kasteleyn orientation given by the "clockwise odd rule" 
[lol [13} around a plaquette. In the thermodynamic limit 
the eigenvalues of K can be used to map out the phase di- 
agram of the model. The phases of the zero temperature 
Kitaev model are therefore in one to one correspondence 
with those of classical dimer models. 

Dimer models divide into two classes: those defined on 
bipartite and non-bipartite lattices. The hexagonal lat- 
tice is bipartite while the Fisher lattice is non-bipartite. 
The physics of the bipartite and non-bipartite lattices are 
quite distinct, so we will discuss them separately. 

We begin with the bipartite case for which one can 
write K in the form 

*-(-5ri) < 5 > 

where the matrix A has real entries and maps between 
the black and white sites. 

Fourier transforming K one obtains a 2d x 2d matrix 
where d is the number of sites of the same colour in a fun- 
damental domain and whose determinant -D(k) = \P\ 2 
where P = Det[A(k)] with k = (0, 4>). 

P, the determinant of A, is then a (Laurent) poly- 
nomial P(e ie ,e 1 ^) in two unimodular complex variables. 
One can move off the unit circle by using some of the cou- 
plings to define general complex variables z and w. We 
will refer the resulting polynomial as the spectral poly- 
nomial of the model. It can generically be written in the 
form 

P(z,w) = ^(-l) a+b+ah p ab z a w b (6) 

a , b 

where a and b are integers and the coefficients p a b are 
positive real numbers determined by the couplings of the 
model. 

The zero locus, 

P(z,w) = 0, (7) 

allows one to immediately identify the gapless phase 
D(k) = 0. This zero locus is called the spectral curve. 




FIG. 1: On the left is the basic tile for the hexagonal lattice, 
showing the activities and our choice of Kasteleyn orientation, with 
the corresponding amoeba on the right. The regions exterior to the 
curves shown corresponds to the Kitaev A phase while the interior 
or amoeba gives the B phase. 



Kenyon, Okounkov, and Sheffield [l4[ and Kenyon and 
Okounkov [15[ proved that the spectral curve |7]), of 
a d x d fundamental domain with generic couplings is 
a Harnack curve and furthermore that every Harnack 
curve arises as the spectral curve of some bipartite dimer 
model. Such curves arose in the early twentieth century 
mathematics literature. They have very special proper- 
ties and can be characterized by what is referred to as 
the amoeba of the curve. 

Given the correspondence of the quadratic form of the 
vortex free sector of the Kitaev model to the Kasteleyn 
matrix of a classical dimer model, as pointed out above, 
cqn. iff]), determines the phase diagram of the zero vor- 
tex sector of the Kitaev model. This phase diagram can 
conveniently be plotted in terms of the amoeba of eqn. 
([7]) as we now describe. 

Put simply, with z = e x+lB and w = e v+l< ^, eqn. (J7J) can 
be solved for the angles 9 and <j> (see [Hj]). There are only 
two solutions (Q(x,y),$(x,y)) and (— &{x, y), — &{x, y)) 
corresponding to (z,w) and (z, w). At the boundary of 
the zero locus these angles become either or n. The 
interior of the domain in the x, y plane bounded by these 
curves is the amoeba. 

In brief the amoeba can be plotted by noting that on 

it 

1 

p,q=0 

The amoeba has bounding curves called ovals. Ovals 
never intersect and can be both compact and non- 
compact. The genus, g, of the spectral curve is equal 
to the number of compact ovals. 

The phase diagram of both classical bipartite dimer 
models and the zero temperature, zero vortex sector of 
Kitaev models then consists of the amoeba and its com- 
plement in the plane. The phase transitions of these 
models occur on the boundary of the amoeba, i.e. on 
the ovals. 

For the lxl fundamental domain with positive cou- 
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FIG. 2: The hexagonal tiling which reduces to the square lattice 
on bond contraction. Dashed lines correspond to bonds that are 
set to zero while all other bonds are allowed to have generic values. 
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FIG. 3: The most general square lattice tiling with two vertices 
of each colour per tile and its amoeba with t z = t w = and 
A = B = 1 and D = 10 in eqn. lUTi 



plings J x , J y and J z the spectral curve is given by 



P(z,w) = 1 - 1/z - 







(9) 



where z = (J y /J x )e and w = (J z / J y )e 1 ^ . We show the 
fundamental tile and its corresponding amoeba in Fig. [TJ 
see [li} for more details on the amoeba in this case. 

The complement of the amoeba consists of both com- 
pact and non-compact regions. In the terminology of 
dimer models, as models of melting crystals, the non- 
compact regions exterior to the bounding ovals consti- 
tute the frozen regions. The amoeba itself is referred to 
as the liquid phase, and the interior of the compact ovals 
as the gaseous phase. 

The amoeba corresponds to the critical surface and 
gives the gapless B phase of the Kitaev model. The 
frozen dimer regions go over to the Kitaev A phase. 
There is however a third C phase corresponding to the 
gaseous phase of the dimer model. 

A gapped phase. — This C phase seems not to have been 
observed in the literature. It is a novel phase in that it 
has a gap which corresponds to a Dirac mass for the 
pair of Majorana Fermions that describe the low energy 
excitations of the model. 

The simplest example of a model that exhibits a C 
phase is the square lattice as shown in Fig. [3l This can 
be obtained, by bond contractions from a hexag- 

onal lattice, with a 3 x 3 fundamental cell, as shown in 
Fig. [2 The Hilbert space of then grows as 2 18rl with 
n the number of copies of the fundamental cell and makes 
direct or Monte Carlo analysis of the phase diagram dif- 
ficult. 

The Kasteleyn matrix takes the form ([5|) and for the 
square lattice with four site fundamental tile, which is 
the simplest model with a compact oval, we have 



A = 



_ ( Jll — JllxS 18 J\1 — Jl2yG 



-i<fj 



(10) 



The determinant of A gives p s i uare which can be pa- 
rameterized in the form 

pSquarc = D _ A , + \ _ B , + _n 

Z W 



where D = 2Acosh{t x ) + IB cosh(t l/ ) and z = e x+ and 
w — e y+l ^ so that z and w can now take any complex 
values. The key point here is that this particular model 
has three phases: Those described in Kitaev's paper [l[ 
as the A and B phases together with a new C phase 
involving a mass gap. This phase arises when t x or t y 
is greater than zero. The phase is bounded by the curve 
P Squaro (e a: , e y ) = which is the compact oval. 

Setting J\\ x = e _t , J22 = e* and all others to unity 
puts the model in its C phase and corresponds to the 
origin of Fig. [3] where t = ln(4 + \/l5). The masses of 
the Majorana modes are mi = 1 — e _t and 7712 = e* — 1. 
For small t the mass gap is given by t and vanishes when 
t = where the C phase disappears. More generally the 
gap can be obtained from the lowest eigenvalue of the 
Kasteleyn matrix, or estimated from the spectral poly- 
nomial whose modulus squared is the determinant of K. 

The general bipartite lattice — all of which can be ob- 
tained by special reductions, involving bond contractions 
ToL HI , of the hexagonal lattice with d x d fundamental 
domain containing d 2 vertices of each colour — is as fol- 
lows: The phase diagram of the zero temperature, zero 
vortex Kitaev model is specified by the amoeba of the 
spectral polynomial which is given as the determinant of 
A as in eqn. (fTTj) . The Kasteleyn matrix, K, then de- 
scribes a set of Majorana Fermions. On the amoeba the 
system is gapless with two of these Fermions being mass- 
less. These in turn become massive as one moves inside 
a compact oval and the system becomes gapped. 

There will in fact be g = (d — l)(d — 2)/2 compact 
ovals. The number, g, of compact ovals is a topological 
invariant and is the genus of the spectral curve. As the 
couplings are varied some of the ovals may contract to 
points. When bond couplings are set to zero, bounding 
phases will disappear. One can remove the corresponding 
bond from the tiling. 

When the system is wrapped on a torus there is an 
additional term in the Hamiltonian 0] and a projection 
operator in the expression for the partition function ([2]). 
However, one can make the following observation: The 
spectrum of the Kasteleyn matrix, K, is discrete and on 
the amoeba the partition function depends on the shape 
parameter, r, (often referred to as the modular param- 



eter) of the torus. As one approaches the non-compact 
boundary of the amoeba the torus becomes degenerate 
and there is a topological transition where the volume 
of the torus goes to zero This can best be seen 

in the context of classical dimers where the finite size 
corrections correspond to the conformal field theory of a 
massless Dirac Fermion. 

Non-bipartite models. — The classical dimer model is 
also useful in gaining an understanding of the Kitacv 
model on a non-bipartite lattice. Again the minimum en- 
ergy configuration will correspond to u std being a Kaste- 
leyn orientation and K(u std ) being a Kasteleyn matrix 
and its critical surface will describe the gapless phase. 

A classic example is the dimer model on the Fisher 
lattice, Fig. [4j whose Kasteleyn matrix is given by 
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(12) 



The Pfafnan of K gives t he p artition function of the two 
dimensional Ising model [lCj . Setting A, B or C to zero 
renders the lattice bipartite and equivalent, by a lattice 
reduction lij, lSj . to the hexagonal lattice. When C = 



P 



= -aA 2 e l 



- bB 2 ^ 



abce 1 ^. 



(13) 



This crossover from bipartite to non-bipartite Ising like 
behaviour was studied in 



17|. 



One can apply this result to the Kitaev model on the 
Fisher lattice, studied in 0. Off the Ising critical surface 
the Kitaev model is gapped and is gapless at the critical 
point. The mass of the lightest Majorana Fermion gives 
the gap. 

Setting A = B = C = 1 and a — 1/vi, b — l/v 2 and 
c = 1/^3, with Vi — tanh(fci), the classical dimer model 
gives the partition function 



^Ising 



s~i yFishcr 
Worm ^ dimer 



, u 2 > "3 



\ 1,1,1) (14) 



for the two dimensional Ising model on the hexagonal 
lattice. The normalization constant C norm is given by 



= 2 



2NM 



sin(fci) sin(fc2) sin(A;3)] 



NM 



(15) 



and Zj^° r r is the dimer partition function on the Fisher 
lattice. 

The hexagonal lattice Ising model (and hence the Ki- 
taev model) is critical for V1V2 + V2V3 + V3V1 = 1; e.g. 
for a = b = c = \/3 with A = B = C = 1, the model is 
gapless with two of the six eigenvalues zero at 9 = <p = n, 
in agreement with 0]. 

Summary. — In this note we established an equivalence 
of the phase diagrams of the vortex free phases of Ki- 
taev models with those of classical dimer models. This 



FIG. 4: The Fisher Lattice. 



allowed us to draw the following conclusions: For bipar- 
tite models the B phase, described by a massless Dirac 
Fermion, corresponds to the amoeba of a classical dimer 
model. Changing the parameters of a sufficiently compli- 
cated model one can cross a phase boundary to a new C 
phase where this Fermion acquires a mass. Adding non- 
bipartite couplings lifts the mass gap degeneracy of the 
two Majorana Fermions. The zero-temperature ground 
state is then a chiral spin liquid When the mass gap 
corresponding to the lightest Majorana Fermion is very 
small the zero vortex sector is well described by a rela- 
tivistic chiral Fermion. 

This work was partly supported by EU-NCG Marie 
Curie Network No. MTRN-CT-2006-031962. 
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